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SUMMARY 

It  has  been  postulated  by  the  author  that  the  burning  rate 
exponent  for  solid  propellant  combustion  can  be  represented  by 

n  ■  log  { 1/Levis  number) 

The  Levis  number,  Le,  is  a  dimensionless  ratio  of  mass 
diffusivity  to  thermal  diffusivity.  The  Levis  number  correlation 
is  used  in  the  paper  to  explain  burning  rate  versus  pressure  be¬ 
havior  for  various  progressive,  plateau,  and  regressive  propellants. 
Qualitative  prediction  by  the  correlation  of  observed  relationships 
of  heat  and  mass  transport,  density,  specific  heat,  and  thermal 
conductivity  is  described.  Discussion  is  also  presented  on  use 
of  the  model  to  predict  a  critical  density  effect  in  deflagration- 
to-detonation  phenomena. 

Using  a  selected  PBAA-ammonium  perchlorate  propellant  as  a 
base,  the  burning  rate  exponents  of  tvo  other  similar  propellants 
vere  calculated.  Calculated  exponent  values,  using  the  Levis  num¬ 
ber  correlation,  vere  vithin  10  to  25  percent  of  measured  exponents. 


TABLE  OF  CONTENTS 


Summary . . . . . 

Table  of  Coutents  . . 

List  of  Figures  . . . . . . . 

Introduction  . . . 

Mathematical  Curves  . . . . . 

A  Hypothetical  Propellant  . . 

Realistic  Levis  Numbers  . 

Cross-Sections  of  Burning  Solid  Propellants 

Qualitative  Prediction  of  Exponent  . 

Quantitative  Prediction  of  Exponent  . 

Deflagration-to-Detonation  . 

Conclusions  . 

References  . . . . 

LIST  OF  TABLES 

1.  Relationship  oi‘  Levis  Number  to  Burning 

Rate  Exponent  . . . 

2.  Propellant  Compositions  . . . 

3.  Comparison  of  Measured  and  Calculated 

Burning  Rate  Exponents  . 

U.  Effect  of  Variables  in  Levis  Number  on 
Deflagration-to-Detonation  . . . 


1 

11 
iii 

1 

2 

J 

3 

5 

6 
9 

12 
15 
IT 


19 

20 

21 

22 


ii 


ITv 


LIST  OF  FIGURES 


1.  Effect  of  Levis  Number  on  r-p  Relationships  . . . 23 

2.  Normal  Rate-Pressure  Relationship  of  Propellants  .............  2b 

3.  Rate-Pressure  Relationship  of  Plateau  Propellants  . . 25 

k.  Rate-Pressure  Relationship  of  Mesa  Propellants  . .  2 6 

.  Hypothetical  Rate-Pressure  Diagram  to  Show 

Different  Propellant  Characteristics . . . . . 27 

6.  Burning  of  a  Solid  Monopropellant  . 28 

7.  Structure  of  Boundary  Region  Associated  with  a 

Burning  Propellant  . . . . 29 

8.  The  Combustion  Region  Adjacent  to  the  Burning  Surface 

of  a  Solid  Propellant  (80*  AP  -  20*  PBAA)  .  30 

9.  Effect  of  Initial  Temperature  on  Rate-Pressure  Relationship  ..  31 

10.  Variation  of  Burning  Rate  over  a  Wide  Pressure  Range  . 32 

11.  Hugoniot  Diagram  . .  33 


THE  PREDICTION  OF  THE  BURNING 


RATE  EXPONENT  OF  SOLID  PROPELLANTS 

BY 

ROBERT  J.  HEASTON 


INTRODUCTION 


The  combustion  of  solid  propellants  involves  a  complex  coupling 
of  transport  phenomena  and  chemical  kinetics.  This  paper  attempts 
to  identify  and  separate  the  effects  caused  by  each  of  these  mechan¬ 
isms.  The  same  model  is  also  used  to  predict  the  characteristics 
and  values  of  the  burning  rate  exponent.  A  unique  approach  to  this 
analysis  is  offered  by  the  following  equation: 


.  loq(VLe)  ,a\ 

r  =  bp^  Cl) 

where:  r  s  regression  rate 

b  s  pre- exponential  function 
p *  pressure 
Le  *  Lewis  number 

./>Ocp 

k 

jo*  density 
D»  mass  dlffusivlty 

specific  heat  at  constant  pressure 
k«  thermal  conductivity 


The  author  first  proposed  equation  1  in  correspondence  to 
Green(l)  and  0rdahl(2).  Some  preliminary  conclusions  concerning 
equation  1  were  published  in  a  review  paper  by  the  author(3). 

Further  analysis  of  equation  1  Indicates  that  it  is  an  Interesting 
framework  for  separating  the  variables  observed  in  solid  propellant 
combustion. 

This  paper  is  divided  into  sections.  Each  section  will  begin 
with  some  prediction  that  may  be  made  with  equation  1.  Immediately 
following  each  prediction  made  by  equation  1  will  be  a  discussion 
of  related  data  from  the  literature.  A  transition  will  then  be  made 
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to  the  next  section.  Over-all  conclusions  of  the  model,  predictions 
and  available  proof  are  given  at  the  end  c-f  the  paner. 


MATHEMATICAL  CURVES 


Independent  of  any  relationship  to  solid  propellant  combustion, 
equation  1  is  a  mathematical  formula  that  can  be  graphed  when  appropri¬ 
ate  numerical  values  are  inserted  for  the  variables.  In  logarithmic 
form  equation  1  is  as  follows: 


log  r  =  log  b  +  (log  1/Le)iog  p 


(2) 


In  a  plot  of  log  r  versus  log  p,  the  value  of  log  1/Le  will  control 
the  slope  or  curvature  of  the  resultant  curves.  Thus,  for  arbitrary 
values  of  r,  p,  and  b,  figure  1  shows  the  curves  for  various  Lewis 
numbers . 

For  Lewis  numbers  approaching  zero,  the  value  of  r  in  figure  1 
is  increasing  infinitely  with  increasing  p.  This  condition  is  de¬ 
pendent  upon  the  fact  that  the  Lewis  number  is  changing  with  p  (or  r 
for  the  purely  mathematical  cese). 

The  widely  known  form  of  equation  1  is  the  de  Saint  Robert  equa¬ 
tion^)  , 


r  =  b  pn  ^ 

Equation  3  is  an  empirical  formula  where  the  exponent  n  is  determined 
from  experimental  measurements  on  actual  propellants.  Different  pro¬ 
pellant  formulations  give  different  regression  rate  versus  pressure 
relationships.  Most  propellants  yield  the  characteristics (l«)  shown 
in  figure  2. 

In  this  case  n  is  usually  greater  than  zero  but  less  than  one. 

A  constant  slope  is  also  maintained  o'er  a  wide  pressure  range  (a  few 
pounds  per  square  inch  to  several  hundred  nounds  per  square  inch). 

Some  nitrocellulose  propellants  with  small  amounts  of  lead  s  >ii- 
tives  nroduce  the  curve(U)  shown  in  figure  3.  As  a  result  of  itF 
characteristic  shane,  the  formulation  producing  the  curve  in  figure  3 
is  called  a  "plateau"  propellant.  A  few  ammonium  perchlorate  composite 
propellants  also  exhibit  this  behavior. 
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Another  unique  relationship  of  regression  rate  and  pressure 
is  exhibited  by  so-called  "mesa"  propellants (4) .  This  figurative 
name  should  be  obvious  from  figure  1*.  A  negative  effect  of  the 
exponent  n  exists  over  a  narrow  pressure  range. 

In  figures  2,  3,  and  1»,  the  values  of  n,  from  equation  3,  are 
preferably  less  than  one.  The  usual  magnitudes,  except  for  negative 
values,  are  Q.U  to  0.9.  A  low  level  is  desirable.  When  n  approaches 
or  exceeds  one,  there  is  a  definite  danger  of  deflagration-to-deto- 
nation  (DTD).  In  DTD,  propellant  regression  increases  almost  infini¬ 
tely  fast  with  pressure  until  a  steady  state  detonation  is  reached. 
Under  these  circumstances  the  exponent  n  is  probably  some  function 
of  pressure. 

It  is  readily  apparent  that  comparison  of  equations  1  and  3 
yields  the  followings 

n=  log  (1/Le)  (4) 

On  a  strictly  numerical  basis,  all  of  the  slopes  depicted  in  figures  2 
through  U  also  are  shown  in  figure  1.  This  relationship  is  the  sub¬ 
ject  of  the  next  section  of  this  paper. 


A  HYPOTHETICAL  PROPELLANT 


A  hypothetical  composite  of  figures  1  through  1*  is  shown  in  fig¬ 
ure  5.  Four  different  zones  are  shown  in  figure  5.  Zone  1  represents 
the  normal  propellant  shown  previously  in  figure  2.  Using  observed 
values  of  n  of  0.1-0. 7,  the  Lewis  number  can  be  calculated  from  equa¬ 
tion  4.  The  resultant  values  are  Le  -  0.8-0. 2.  In  zone  2  of  figure 
5,  Le  «*  1.0  for  n  ■  C.  Similar  determinations  may  be  made  for  zones 
3  and  1*.  Thus,  all  the  values  of  the  Lewis  number  presented  in  figure 
l.also  occur  in  figure  5.  Every  magnitude  of  n  reported  in  figures 
2-1*  is  repeated  in  figure  5.  With  the  exception  of  zone  1*,  figure  5 
is  very  much  like  figure  L.  Instead  of  changing  to  a  normal  propel¬ 
lant  slope  in  zone  1*  in  figure  5f  a  deflagration-to-detonation  cha¬ 
racteristic  is  substituted.  The  mesa  of  zone  2  is  also  sore  idealistic 
with  a  horizontal  behavior.  It  is  quite  possible,  then,  that  figure  5 
is  not  so  hypothetical  after  all. 


REALISTIC  LEWIS  NUMBERS 


The  values  of  the  Lewis  numbers  shown  in  figures  1  and  5  must 
have  some  physical  significance  or  else  the  hypotheses  made  so  far 
would  be  meaningless.  The  Lewis  number  is  a  dimensionless  ratio  that 
is  often  interpreted  as  the  Prandtl  number  divided  by  the  Schmidt  number. 
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3H 

ii 

CL 

Prandil 

number 

Sc  =  JL. 
/JO 

Schmidt 

number 

(5) 

(6) 


In  these  dimensionless 
previously,  except  for 
(6), 


Le  = 


Pr 

Sc 


numbers  the  terms  are  the  same  as  defined 
M  which  is  viscosity.  Dividing  (5)  by 

(7) 


Using  values  of  Pr  and  Sc  in  Foust  «t  al ( 5 ) .  the  values  of  Le  shown 
in  table  1  are  matched  with  a  description  of  real  gases.  Real  gases 
generally  have  Lewis  numbers  nearly  equal  to  unity  or  slightly  less. 
The  Lewis  numbers  of  liquids  or  gases  at  very  high  pressures  are 
usually  quite  low.  For  the  ideal  case  using  a  simple  model  gas,  the 
transport  diffusivities  are  equal  (5).  Since,  alpha,  the  thermal 
diffusivity  is 


a  = 


k 

jOCp 


This  means  then,  that 


Le  = 


D. 

a 

Mass  diffusivity 
Thermal  diffusivity 


(8) 


(9) 


As  stated  above,  the  mass  diffusivity  and  the  thermal  diffusivity 
are  eoual  for  sample,  ideal  gases.  Consequently,  Le  ■  1  for  an 
ideal  gas. 

Unfortunately  Levis  numbers  are  difficult  to  measure  in  real 
gas  systems.  It  becomes  even  more  difficult,  if  not  impossible, 
to  define  Lewis  numbers  for  a  complex  system,  especially  a  burning 
solid  propellant. 
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CROSS-SECTIONS  OF  BURNING  SOLID  PROPELLANTS 


In  order  to  be  a  useful  tool,  the  Levis  number  must  be  defined 
for  some  part  of  the  combustion  zone  that  exists  in  burning  solid 
propellants.  The  problem  is  even  more  difficult  because  th*re  is 
some  variation  in  the  view-points  of  vhat  actually  occurs  as  a  pro¬ 
pellent  burns.  £  lid  propellants  are  usually  a  heterogeneous  mixture 
of  fuel,  binder,  and  oxidizer.  In  an  organic  composite  propellant, 
these  ingredients  are  usually  powdered  aluminum,  polymer,  and  ammonium 
perchlorate,  respectively.  A  doubl^-baae  propellant  might  contain 
powdered  alumin  .n  in  a  plasticized  nits,  -t4-,  ...lose  formulation.  Var¬ 
ious  additives  are  used  t~  obtain  difium.'v  ret  i -pres sure  relation¬ 
ships.  I’nder  certain  coriu  tions  one  might  also  argue  whether  tho 
oxidizer  regresses  a  little  faster  than  the  binder,  or  that  a  molten 
xayer  exists  on  the  burning  propellant  surface.  Nevertheless,  some 
general  agreement  has  bean  reached  on  vhat  occurs.  Three  separate 
crcs&  sections  perpendicular  to  the  burning  propellant  surface  are 
shown  in  figures  6,  7,  and  8. 

Figure  6  presents  nomenclature,  such  as  "fizz  zone"  and  "foam 
zone,"  that  has  been  observed  to  behave  in  the  terms  described(U) . 

A  mechanistic  cross  section  is  shown  in  figure  7.  This  model  by 
Cantrell  et  al.  (6)  also  includes  combustion  instability  effects. 

A  more  physical  interpretation  (7)  of  the  burning  surface  of  a  solid 
propellant  is  shown  in  figure  8. 

As  reported  by  Vantock{8),  the  Soviets  also  differentiate  stages 
in  the  solid  propellant  combustion  process.  These  stages  are:  (a) 
preheating  of  the  propellant;  (b)  solid  phase  reactions;  (c)  melting; 
(d)  liquid  phase  reactions;  (e)  evaporation;  (f)  sublimation  or  dis¬ 
persion;  and  (g)  homogeneous  or  heterogeneous  reaction  in  the  gae 
phase.  Regardless  of  which  model  (whether  figure  b,  7,  or  8)  is 
chosen  it  should  be  obvious  that,  a  complex  process  occurs  across  the 
interface  between  hot  combustion  gaais  and  the  unreacted  solid  pro¬ 
pellant.  For  the  time  being  it  must  be  assumed  that  a  real  Lewis 
number,  or  an  "effective”  Lewis  number,  can  be  defined  for  the  over-all 
combustion  region. 

It  is  significant  to  note  that,  rather  than  acknowledge  the 
existence  of  such  a  complex  function  as  a  so-called  "effective"  Lewis 
number,  the  normal  practice  is  to  assume  Le  »  1.0.  However,  there 
have  been  evaluations  of  situations  when  Le  t  1.0.  One  of  the  be3t 
of  these  analyses  is  Lee8(9).  In  boundary  type  How  systems  with 
chemical  reactions,  heat  energy  is  transported  bv  heat  conduction 
and  the  diffusion  of  species  carry'. ,g  chemical  enthalny .  A  heat,  flux 
vector,  4,  can  be  defined  to  represent  thin  over-all  heat  transport. 


(10) 


Thus , 


-Jslfh  -2h)|Ki  .^%?PEh||Kl 
CpLlay  dyl  k  9y 


CONDUCTION 


DIFFUSION 


where:  hj  =^cpjdT  ♦  h* 

rfj  s  heat  of  formation  of  f'th  species 
h  ■  complete  static  enthalpy 
=  EK(h| 

K|  *  mass  fraction  of  i'th  species 
D,2a  binary  mass  dlffuslvlty 

Cp  a  £h|Cp| 

T  =  temperature 


When  the  Levis  number  equals  unity  in  equation  10,  the  heat  flux 
vector  becomes 


4=  - 


_k  3h 
cP  dy 


Equation  11  is  similar  tc  ordinary  heat  conduction  for  a  non-reacting 
pure  gas.  This  is  usually  a  much  easier  equation  to  sol va  than  equa¬ 
tion  10.  Even  though  the  Levis  number  may  be  eliminated  from  calcu¬ 
lations  by  assuming  it  equal  to  unity,  there  is  no  overlooking  the 
fact  that  the  Levis  number  haB  significance  in  physical  phenomena. 

A  qualitative  interpretation  of  the  Levis  number  and  equation  1  vill 
indicate  some  of  the  phenomena  involved  in  r.olid  propellant  combustion. 


QUALITATIVE  PREDICTIOH  07  EXPOHSKT 

Baaed  upon  equation  1,  and  experimental  evidence,  there  are  four 
major  parameters  that  '-.^tribute  to  differences  in  solid  propellant 
burning  rates  at  constant  pressure.  These  arei 

4.  Chemical  Composition 

b.  Transport  Processes 

c.  Chemical  Kinetics 

d.  Physical  Ovaetry 
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Parameters  a,  b,  and  c  can  be  partly  related  through  use  of  figure  5. 
Assume  a  propellant  is  burning  with  the  conditions  and  properties  a3 
shown  in  zone  1  of  figure  5.  If  this  hypothetical  propellant  is  like 
a  plateau  or  mesa  propellant,  its  rate-oressure  curve  will  change  in 
the  neighborhood  of  some  pressure  "x".  What  happenn  at  x  (which  does 
not  necessarily  have  to  be  a  sharp  transition)  ?  Supposedly,  the  ki¬ 
netic  order  of  the  system  changes  between  zones  1  and  2.  In  terms 
of  the  Lewis  number  the  transition  involves  a  coupling  between  trans¬ 
port  properties  and  kinetic  phenomena.  In  zone  1,  it  is  hypothesized 
that  reaction  A  — ►  B  predominates  in  the  combustion  process  near 
the  propellant  surface.  Under  these  conditions  the  Lewis  number  will 
remain  constant  and  largely  dependent  upon  the  species  B.  As  pressure 
increases  to  zone  2,  then,  reaction  A  — C  becomes  significant.  A 
new  species,  C,  now  controls  the  flame  properties  and  a  new  Lewis  num¬ 
ber  is  established.  If  the  reaction  does  not  change,  the  same  Lewis 
number  should  predominate  until  increased  pressure  starts  affecting 
the  transport  processes.  Spectroscopic  observations  of  the  com¬ 
bustion  zones  of  different  types  of  propellants  burning  at  different 
pressures  could  possibly  prove  or  disprove  this  mechanism.  Some 
data  along  this  line  have  been  reported  by  Tajima  et  al.(lO).  In 
these  studies,  spectroscopic  data  indicate  that  the  same  chemistry 
is  taking  place  in  the  surface  reaction  zone  of  double-base  pro¬ 
pellants  regardless  of  the  pressure.  No  change  in  burning  rate  ex¬ 
ponent  was  reported  in  these  studies. 

In  considering  figures  2  through  4,  it  should  be  obvious  in  the 
light  of  present  experience,  that  certain  propellant  compositions 
have  inherent  burning  characteristics.  The  problem,  then,  becomes 
cne  of  how  to  modify  these  characteristics. 


Combining  equations  4  and  9  the  following  relation,  'ip  of  burning 
rate  exponent  and  transport  processes  is  obtained: 


n  = 


<t> 


Therm  a  I  _  d  i  t  j  u  s  i  v  it  y 
Mass  diffusivity 


(12) 


Consequently  the  burning  rate  exponent  will  be  increased  by  increases 
in  the  thermal  diffusivity  if  the  mass  diffusivity  remains  unchaiged. 
Increased  mass  transport  will  lower  the  exponent  if  no  changes  in 
thern_.l  transport  occur.  However,  it  is  difficult  (if  not  impos  ible) 
to  change  the  transport  processes  without  changing  the  kinetic  m  cha- 
nism.  It  has  been  noted(4)  that,  at  rocket  pressures  (below  20C  '  psi), 
n  for  the  same  propellant  is  generally  lower  than  at  gun  pressures 
(10,000  to  50,000  nsi).  This  is  probably  due  to  the  fact  that  higher 
pressures  tend  to  increase  thermal  diffusivity  and  retard  mass  dif¬ 
fusivity. 
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As  reported  by  Vantoch(lO),  it  has  been  predicted  and  observed 
that  combustion  is  stabilized  with  an  increasing  amount  of  heat 
released  in  the  condensed  phase.  This  vould  lead  to  a  lowered 
thermal  diffusivity  and  a  correspondingly  lower  or  more  stable 
exponent.  Vantoch(lO)  also  reports  that  the  effect  of  a  catalyst 
in  some  cases  is  to  cause  a  high  surface  temperature  through  heat 
release  in  the  condensed  phase.  Heat  transfer  from  the  gas  to  the 
condensed  phase  would  have  to  be  small  under  these  circumstances. 

Another  approach  would  be  to  cause  an  apparently  artificial 
increase  in  the  thermal  diffusivity  through  external  radiation  or 
an  additive  that  considerably  increases  the  flame  radiation.  This 
could  be  done  without  appreciably  disturbing  the  chemical  steps  in 
the  reaction  mechanism.  A  small  amount  of  heavy  metals  might  also 
disturb  the  mass  diffusivity  sufficiently  to  lower  the  exponent. 
Negative  exponents  are  often  observed  when  condensed  residues  from 
heavy  metals  or  other  additives  are  formed  during  combustion.  When 
a  condensed  residue  forms,  the  regression  rate  decreases  with  in¬ 
creasing  pressure  due  to  lowered  heat  transfer  (10). 

The  orientation  of  metal  fibers  in  a  prooellant  could  also 
cause  increased  thermal  diffusivity,  resulting  in  higher  exponents. 
However,  this  effect  also  would  depend  upon  the  kinetics  of  the 
metal  combustion  process.  If  the  metal  fibers  react  at  a  rate  less 
than  the  unmetallized  propellant,  the  effect  of  the  increased  ther¬ 
mal  diffusivity  would  be  the  same  as  a  higher  ambient  temperature 
in  the  propellant.  In  this  case,  the  over-all  burning  rate  is  in¬ 
creased  with  little  or  no  change  in  burning  rate  exponent. 

In  equations  1  and  3#  the  quantity,  "b” ,  is  the  variable  which 
depends  primarily  on  the  chemical  kinetics.  This  conclusion  is 
most  readily  apparent  from  the  effect  of  initial  temperature  on 
propellant  burning  rates.  Figure  9  illustrates  this  point  (U). 
However,  with  equation  1  it  becomes  more  evident  that  b  is  dependent 
on  kinetics.  The  same  effect  as  produced  by  variation  in  initial 
propellant  temperature  is  also  accomplished  by  ihe  addition  of  ca¬ 
talysts  to  solid  propellants.  The  exponent  might  be  changed  by  adding 
a  sufficiently  large  quantity  of  catalyst  that  the  transport  pro¬ 
cesses  are  affected.  Catalysts  which  ehange  the  intermediate  reaction 
steps  might  also  cause  a  variation  in  exponent.  After  all,  any  change 
in  composition  or  the  geometry  of  constituents  means  a  new  propellant. 

The  effects  of  oxidizer  particle  size,  propellant  homogeniety, 
and  physical  geometry  have  been  noted  by  Hall  and  Bastress(ll) . 

The  results  in  general  show  the  same  variation  as  the  effect  of 
initial  temperature  or  catalysts. 
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Throughout  this  section,  most  of  the  comparison  of  propellant 
properties  have  been  made  at  a  constant  pressure.  However,  pressure 
can  be  a  significant  variable,  too.  As  mentioned  before,  the  changes 
that  occur  in  figure  5  at  points  x,  y,  and  i  are  a  function  of  pres¬ 
sure  that  are  related  to  the  intrinsic  properties  of  the  propellant. 
There  is  another  significance  to  pressure  that  can  best  be  identified 
by  another  hypothetical  model  propellant .  These  properties  are  shown 
in  figure  10.  Any  of  the  typical  propellant  performance  curves  siiown 
in  figures  2,  3,  and  ^  could  have  been  depicted  in  the  "normal  pro¬ 
pellant  operations"  region  of  figure  10.  The  propellant,  represented 
by  figure  2,  is  a  more  common  example.  Pure  ammonium  perchlorate 
and  most  propellants  made  with  this  oxidizer  show  a  nearly  linear 
behavior  over  an  extremely  wide  pressure  range.  In  fact,  rate- 
pressure  measurements  on  a  polybutadiene-acrylic  acid  composite 
propellant  have  been  made  by  Cole(l2)  up  to  60,000  psi  in  a  self- 
pressurizing  closed-bomb  system.  These  studies  did  not  exhibit 
a  deflagration-to-detonation  transition. 

At  very  low  pressures,  the  combustion  zone  broadens.  Diffusion 
processes  become  controlling  to  the  point  where  an  extinction  pres¬ 
sure  is  reached.  Equation  1  conforms  with  such  an  approach. 

Under  certain  conditions  when  confined  in  a  bomb  or  subjected 
to  an  energy  transferring  impulse,  the  pressure  rises  at  a  very 
high  rate.  A  transition  then  occurs  from  normal  deflagration  to  a 
steady  state  detonation.  If  not  suitably  confined,  the  pressure 
rise  may  only  result  in  a  catastrophic  explosion.  Typical  shock 
pressures  for  propellants  may  be  60  kilobars  with  "regression  rates" 
of  5000  meters  per  second.  However,  before  deflagration-to-detona- 
tion  processes  can  be  described  in  more  detail,  it  is  necessary  to 
discuss  the  specific  effect  of  variables  on  the  burning  rate  exponent. 


QUANTITATIVE  PREDICTION  OF  EXPONENT 

It  would  be  desirable  to  be  able  to  predict  the  magnitude  of  the 
burning  rate  exponent  of  a  propellant  from  its  physical  properties. 

If  this  can  not  be  done,  it  might  also  be  just  as  desirable  to  predict 
the  properties  of  one  propellant  from  the  properties  of  another.  The 
following  analysis  considers  the  latter  point. 
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Using  equation  1  and  comparing  propellant  1  against  propel¬ 
lant  2, 


n  = 

,  loqO/Le*) 

bi  Pi 

(13) 

ti 

log(1/Le2) 

b2  p2 

(14) 

For  relative 

rates  at  equal  pressures. 

r* 

-U  - 

Hi 

b1  logd/Le^  -  log(1/Le2) 

(15) 

rvP  Lr> 

II 

b^  log  (Le2/  Le^) 

(16) 

Make  use  of  the  following 

ax  = 

In  a 

e 

(17) 

In  b 

a 

In  b  (in  a) 

=  e 

(18) 

This  is  the  same  as 

,  In  a  (In  a)  In  b 

b  =  e 

Then  t 


a 


I  n  b 


or. 


glog  b  _ 


k'oq  a 


(19) 

(20) 

(21) 
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equation  (16)  so  thatj 


The  identity  above  can  be  used  in 


Since  y 


Assume 


n 

P2  ’ 

b, 

b2 

Le2' 

iLe'i 

!og  p 

Le  = 

/>Dcp 

k 

_  /  ^°2  D2Cp  \ 

n  . 

bi 

k2 

r2  ' 

b2 

fi i  D,cp,  1 

^  1 

y°i 

•  *2  .  < 

Pi 


cp2  ,  Di  -  D2 


ill 

k2 

\  2i 


log  p 


(22) 

(23) 

(24) 


(25) 


Of...  7  \  ine  above  anal' 
"ea!®  in  boning  rate  from  r,  to  i 
elative  changes  in  variables  occuj 


cay  be  summarised.  An  in. 
results  if  the  following 


ra>f?  If:  bj>  b, ,  other  variables  equal 
ka  >k,.  other  variables  equal 
>  other  variables  equal 
other  variables  equal 
S<  S'  other  variables  equal 
^>0^,  other  variables  equal 
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Data  from  a  study  by  Sabadell  and  Wenograd(l3)  were  used  to 
make  calculations  related  to  equations  22  to  25.  In  the  studies 
by  Sabadell  and  Wenograd  the  surface  temperatures  and  thermal 
conductivities  of  a  series  of  PBAA-amnonium  perchlorate  propel¬ 
lants  were  measured.  The  composition  of  the  propellants  of 
interest  are  shown,  in  Table  2.  Composition  70F  was  used  as  a 
base  to  calculate  a  Lewis  number.  The  following  values  were  cal¬ 
culated  or  assumed  as  order-of-magnitude  quantities: 


fi  »  1.8  grams/cm^ 

cp  =  1.0  cal/°C 

log(i-)  =  °‘255 

Le 

k  *  5.3  x  10-4  cai/°C-cm-oec 


p  *  15  psig 

D  1.652  x  10-4  cm2/sec 


Other  measured  values  of  thermal  conductivity  were  substituted 
into  log  (l/Le)  to  determine  the  calculated  exponents  in  Table  3. 

The  calculated  burning  rate  exponents  do  not  agree  too  close¬ 
ly  with  the  measured  values.  However,  the  calculated  quantities 
bracket  the  measured  exponent  for  the  two  different  thermal  con¬ 
ductivities  used.  It  is  significant  to  note  how  the  changes  in 
thermal  conductivity  change  the  calculated  exponents.  Some  other 
properties  such  as  density  or  diffUBivity  must  also  have  changed 
with  the  thermal  conductivity.  Otherwise,  the  exponent  of  the 
individual  PBAA  propellants  would  not  remain  constant  over  a  pres¬ 
sure  range.  In  any  case,  no  other  example  was  discovered  in  this 
review  where  burning  rate  exponents  of  one  propellant  were  pre¬ 
dicted  from  properties  of  another  propellant  formulation.  This 
example  provides  a  crude  but  quantitative  demonstration  of  equa¬ 
tion  1. 


DEFLAGRATION  -TO-DETOMTIOtl 


There  are  circumstances  under  which  a  burning,  or  deflagrating 
propellant  begins  to  regress  at  a  higher  and  higher  rate.  Zone  It 
in  figure  5  represents  this  condition.  Predictions  of  the  mecha¬ 
nisms  of  deflagration-to-dotonation  (abbreviated  as  DTD)  can  be 
made  with  the  aid  of  equation  1.  In  order  to  do  this  it  is  necessary 


12 


first  to  refer  to  figure  11.  The  pressure  versus  inverse  density 
relationship  shown  in  figure  H  is  known  as  a  Hugoniot  plot.  Two 
regions  exist  in  this  plot  with  a  discontinuity  between  the  two 
curves.  As  noted,  one  curve  describe?  conditions  for  detonation 
and  the  other  designates  deflagration. 

In  order  to  correctly  predict  the  influence  of  the  Lewis  num¬ 
ber  on  DTD,  it  is  necessary  to  refer  to  both  figure  5  and  figure  11. 
Consider  the  deflagration  curve  in  figure  11.  As  the  density  gets 
smaller,  the  deflagration  curve  approaches  the  abscissa.  Apparent¬ 
ly,  the  density  could  decrease  to  zero  and  the  deflagration  curve 
could  be  extended  asymtotically.  However,  according  to  equations 
1  and  U,  as  the  density  approaches  zero,  the  burning  rate  exponent 
approaches  an  infinite  magnitude.  This  is  contrary  to  figure  5 
as  well  as  experience.  The  burning  rate  exponent  must  remain  less 
than  unity  if  deflagration  is  to  be  sustained.  A  value  of  n  greater 
than  unity  would  lead  to  a  detonation.  The  transition  of  n  to 
greater  than  unity  occurs  in  equation  k  for  a  Lewis  number  of  about 
0.1.  Thus,  in  figure  11  a  critical  density  is  predicted  where  the 
deflagration  curve  involves  a  transition  to  detonation.  Observations 
of  this  effect  have  been  recorded  by  Gordon(l5).  It  has  also  been 
observed  that  decreasing  the  bulk  density  of  a  propellant,  such  as 
shredding  it,  will  cause  an  otherwise  non-sensitive  material  to  de¬ 
tonate.  If  the  Lewis  number  is  slightly  less  than  unity,  a  small 
decrease  in  density  could  cause  tne  Lewis  number  to  exceed  unity. 

This  would  cause  n  as  defined  by  equation  4  to  yield  a  runaway  re¬ 
action.  It  is  also  possible  that  the  same  process  of  shredding, 
or  high  porosity,  could  lead  to  higher  localized  thermal  diffusivity, 
such  as  a  hot  spot  (l6).  According  to  equation  9,  this  would  also 
lead  to  a  transition  to  detonation.  Thus,  both  the  effects  of  poro¬ 
sity  and  the  potential  existence  of  hot  spots  are  predicted  by 
equation  1. 

In  equation  1,  when  the  Lewis  number  is  unity,  the  regression 
rate  is  apparently  not  a  function  of  pressure.  The  value  of  r  is 
dependant  on  the  kinetic  function  b  only.  Consequently,  in  figure  11, 
the  discontinuity  between  the  detoiation  and  deflagration  curves  pro¬ 
bably  occurs  for  Le  »  1.0. 

The  various  values  of  Lewis  r, umber  and  burning  rate  exponent 
that  exist  for  figure  11  are  tabulated  in  Table  i*.  Zone  3  in  figure 
5  is  also  indicated  under  the  deflagration  curve.  This  is  an  unusual 
case  where  diffusion-related  phenonena  predominate.  The  case  when 
Le  ■  1.0  is  explained  by  the  following  postulation.  Assume  that 
transport  processes  and  chemical  kinetics  can  be  competing  mechanisms, 
even  in  a  detonation.  In  addition,  there  is  always  a  coupling  oetween 
these  mechanisms.  When  a  shock  is  passed  through  an  explosive,  the 
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shock  effects  both  the  transport  processes  and  kinetics.  Suppose 
that  the  Lewis  number  is  unity,  and  is  unchanged  by  a  shock.  The 
processes  contributing  to  reaction  kinetics  are  then  controlling 
when  a  shock  passes  through  the  propellant.  The  shock  must  be 
equal  to  or  greater  than  the  energy  barrier  required  to  initiate 
the  detonation.  If  the  shock  is  not  intensive  enough,  there  will 
be  no  detonation  because  the  Lewis  number  is  unity.  With  a  suffi*- 
ciently  intensive  shock,  a  steady  state  detonation  is  established. 
This  is  the  general  mechanism  of  a  high  order  detonation.  The  de¬ 
tonation  model  supported  by  Eyring  and  coworkers,  as  quoted  by 
Boyer(l7),  follows  this  latter  model.  According  to  these  investi¬ 
gators  the  high  regression  rates  of  a  detonating  grain  in  the  grain 
burning  process  effectively  prevents  conduction  of  heat  into  the 
unburned  solid. 

A  true  deflagration-to-detonation  event  requires  the  Lewis 
number  properties  given  in  Table  U.  There  are  also  circumstances 
where  an  explosive  receives  a  shock  that  is  insufficient  for  a  high 
order  detonation.  However,  the  shock  influences  the  Lewis  number 
such  that  ii  exceeds  unity.  The  reaction  grows  at  a  rapid  pace  be¬ 
cause  of  the  high  pressure  inherent  in  the  initial  shock.  Rapidly 
the  explosive  reaches  a  steady-state  detonation.  But,  this  process 
requires  several  more  microseconds  than  a  high  order  detonation. 

This  slower  process  where  transport  processes  are  significant  results 
in  the  so-called  low  order  detonation  phenomena. 

The  Lewis  number  can  be  unity  for  either  a  detonating  or  de¬ 
flagrating  system.  As  noted  above,  if  Le  »  1.0,  a  sufficient  shock 
must  be  used  to  achieve  a  detonation.  After  a  steady-state  detona¬ 
tion  is  achieved,  the  Lewi3  number  always  adjusts  to  unity.  Other¬ 
wise,  the  reaction  rate  would  increase  indefinitely.  It  might  be 
hypothesized  that  the  mechanisms  for  all  transport-related  phenomena 
become  equivalent  in  detonation  or  supersonic  flow  systems.  For  a 
deflagration,  the  Lewis  number  can  be  any  value  above  about  0.1. 

Some  Justification  for  the  above  mechanism  may  be  provided  by 
combustion  instability  studies.  As  noted  by  Vantoch(lO),  combustion 
instability  may  precede  deflagration-to-detonation  transition.  A 
Dulse  created  through  combustion  instability  may  be  caused  to  grow. 
According  to  Zucrow(l8),  any  effect  which  contributes  to  an  increase 
in  3tatic  temperature  of  an  element  of  gas,  such  as  chemical  heat 
addition,  will  help  to  amplify  a  pulse.  This  could  also  cause  an 
increase  in  thermal  diffusivity  sufficient  for  a  detonation  to  occur. 

Although  not  defined  in  this  paper,  the  existence  of  a  critical 
diameter  for  propellants  and  explosives  is  also  probably  a  function 
of  the  Lewis  number  and  transport  processes. 


CONCLUSIONS 


It  has  been  postulated  by  the  author  that  the  burning  rate 
exponent  for  solid  propellant  combustion  can  be  represented  by 

n  a  log  (l/Lewia  number) 

The  Levis  number,  Le,  is  the  ratio  of  mass  diffusivity  to  thermal 
diffusivity.  Although  this  new  burning  rate  exponent  function  of  Le 
has  not  bsen  rigorously  derived,  evidence  has  been  collected  to  Justi¬ 
fy  some  aspects  of  its  validity.  Such  a  function  must  have  more  than 
Just  a  fortuitous  significance  in  view  of  the  following: 

1.  Measured  values  of  "n”  substituted  in  the  correlation 
produce  reasonable  order-of-magnitude  numerical  values  for  Le. 

2.  The  Le  function  offers  a  single  parameter  as  a  unified 
explanation  for  normal,  mesa,  plateau  and  negative  burning,  and  dei- 
fiagratiou-to-detonation  behavior  of  solid  propellants.  This  is  the 
first  correlation  to  cover  such  a  broad  regime. 

3.  In  comparing  the  burning  rate  properties  of  two  different 
propellants,  the  Le  correlation  predicts  the  observed  relationships 

of  density,  specific  heat,  and  thermal  conductivity. 

U .  Calculations  on  a  selected  propellant  system  gave  an 
order  of  magnitude  prediction  of  the  burning  rate  exponents  of  tvo 
propellants  as  related  to  a  third  propellant. 

5.  The  Le  correlation  complies  with  observed  qualitative 
effects  of  variations  in  heat  and  mass  transport. 

6.  The  Le  correlation  predicts  a  critical  density  effect 
in  detonation  phenomena. 

7.  The  burning  rates  of  solid  propellants  are  expressed 

in  the  Le  correlation  as  a  function  of  heat  transfer  divided  by  mass 
transfer.  It  does  not  make  any  difference  whether  the  numerator  ap¬ 
proaches  infinity  or  the  denominator  becomes  negligibly  small,  the 
burning  rate  increases  very  sharply.  This  essentially  shows  that 
deflagration-to-detonation  analyses  c«n  use  either  mass-diffusivity 
or  heat -resistance  controlled  grain  burning  models. 

8.  A  coupling  effect  between  transoort  properties  and  kinetic 
phenomena  is  exhibited  by  the  Le  correlation.  As  long  as  the  reaction 
A  — B  predominates ,  the  value  of  Le  will  remain  constant  and  de¬ 
pendent  upon  the  3pecie3  "B".  As  pressure  increases,  then,  the  reac¬ 
tion  A  — ►  C  may  become  important.  In  this  case,  a  new  Le  is 
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determined  by  the  speeias  "C".  Such  a  pressure  effect  on  the  con¬ 
trolling  reaction  in  solid  propellant  combustion  has  been  observed, 

9.  The  Le  correlation  allows  a  very  vide  pressure  depend¬ 
ency  which  has  been  a  weakness  in  other  burning  rate  models. 

Even  though  all  of  the  above  may  be  true,  it  is  still  necessary 
to  observe  that  the  "effective  Le"  involved  in  solid  propellant  com¬ 
bustion  may  be  extremely  difficult  or  impossible  to  define,  never¬ 
theless,  the  Le  correlation defined  may  be  an  effective  tool  fir 
analyzing  the  proper  relationships  of  variables  in  solid  propellant 
combustion,  or  heterogeneous  reactions  in  general. 


QjoIUa^  v . 

ROBERT  J.  HE^STOH 
Chief,  Chemistry  Branch 
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TABLE  1 


RELATIONSHIP  OF  LEWIS  NUMBER 
TO  BURNING  RATE  EXPONENT 


Measured 

Exponent 

Rate-Pressure 

Curves 

Calculated 
Lewis  No. 

Description 
of  Gas 

0.1  -  0.7 

Normal  or 
Progressive 

0.8  -  0.2 

Real 

0 

Plateau 
or  Mesa 

1.0 

Ideal 

<0 

Negative  or 
Regressive 

>1 

Complex 

CO 

Deflagration 
to  Detonation 

— ►  0 

Dense 
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TABLE  2 


PROPELLANT  COMPOSITIONS 


Components 

-  Percent  in  Weight 

Propellant 
Dea Ignat Ion 

PBAA 

EPON 

828 

AP  Fine 

9  Micron 

AP 

Coarse 

Ground 

Copper  Ferric 

Chromite  Oxide 

70F 

25.7 

4.3 

70 

...  ... 

70F-1C 

25.7 

4.3 

70 

— 

1 

70F-1F0 

25.7 

4.3 

70 

— 

1 

75  AM 

21.4 

1 

3.6 

22.5 

52.5 

’ -  - 
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TABLE  3 


COMPARISON  OF  MEASURED  AND  CALCULATED  BURNING  PATE  EXPORTS 


i 


Propellant 

Designation 

Pressure 

P«ig 

Thermal 

Conductivity 

cal/°C-cm-sec 

Measured 

n 

Calculated 

n 

70F 

15 

5.3  x  10“* 

0.255 

Bast 

100 

8.4  x  10“* 

0.255 

Base 

TOF-1C 

15 

10.9  x  10"* 

0.381 

0.563 

100 

18.6  x  10-* 

0.381 

0.364 

70F-1F0 

15 

10.6  x  10“* 

0.410 

0.551 

100 

15.3  x  10-* 

0.410 

0.352 

75  BM 

15 

5.5  x  10“* 

0.291 

0.266 
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TABLE  U 


EFFECT  OF  VARIABLES  IN  LEWIS  NUMBER  ON  DEFLAGRATION  TO  DETONATION 


Variable 

Detonation 

Uncertain 

Deflagration 

Le 

0  to  0.1 

1.0 

0.1  to  <  lj  >  1.0 

n 

>1.0 

0 

>0  to  <1;  <0 

fi 

Low 

— 

High 

D 

Low 

— ■ 

High 

CP 

Low 

High 

k 

High 

— 

Low 

at 

High 

Low 
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LOG  p 


Fig.  1  Effect  of  Levis  Number  on  r-p  Relationships 
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LOG  PRESSURE 


i 

i 

t 

! 

[  Fig,  2  Normal  Rate-Pressure  Relati r.nwMp  of  Propellants 


I  2U 


l 


x 


LOG  PRESSURE 


fig,  U  Rata-Praaaura  Ralatlonahlp  of  Maaa  Propallanta 
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RATE 


Fig.  5  Hypothetical  Rate-Pressure  Diagram  to  Show 
Different  Propellant  Characteristics 
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A  COLD  PROPELLANT 

B  HEATED  PROPELLANT 

C  PYROLYSIS  ZONE  (FOAM  ZONE) 

D  PRIMARY  GASEOUS  PYROLYTIC 
PRODUCTS 

E  GAS  HEATING  ZONE  (FIZZ  ZONE) 

F  GAS  REACTION  ZONE  (FLAME  ZONE) 
G  BURNED  GAS 


Fig.  6  Burning  of  *  Solid  Monopropollant 
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BURNING 

ZONE 


PRODUCT 

GAS 


Fig.  T  Structure  of  Boundary  Region  Associated  with  a 
Burning  Propellant 
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PROPELLANT 


£ 


LOG  PRESSURE 


Fig,  10  Variation  of  Burning  Rate  over  a  Wide  Pressure  Range 
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SPECIFIC  VOLUME  C/)'1) 


Fia.  IX  Hugoniot  Diagram 
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